Drift causes anomalous exponents in growth processes 



Gunnar Pruessneqj 

Department of Mathematics, Imperial College London, 180 Queen's Gate, London SW7 2BZ, UK 

(Dated: 23 October, 2003) 

The effect of a drift term in the presence of fixed boundaries is studied for the one-dimensional 
Edwards- Wilkinson equation, to reveal a general mechanism that causes a change of exponents for 
a very broad class of growth processes. This mechanism represents a relevant perturbation and 
therefore is important for the interpretation of experimental and numerical results. In effect, the 
mechanism leads to the roughness exponent assuming the same value as the growth exponent. In 
the case of the Edwards- Wilkinson equation this implies exponents deviating from those expected 
by dimensional analysis. 
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The Edwards- Wilkinson (EW) equation , as it is dis- 
cussed below, is probably the best-studied equation de- 
scribing surface growth processes. Due to its linearity 
it is solvable by standard methods and has been stud- 
ied analytically as well as numerically in great detail 
0- S 13 • The equation is very well-behaved, so that 
the outcome of these investigations are usually quite pre- 
dictable. There is no reason to suspect that well-accepted 
methods, such as dimensional analysis and coarse grain- 
ing, produce wrong results, even if applied to the EW 
equation with an extra drift term (EWd), which still rep- 
resents a linear problem. 

However, it is shown below that, depending on the 
boundary conditions, such a drift term changes the ex- 
ponents dramatically to anomalous values, which appar- 
ently have been missed in the literature. The drift in 
conjunction with the boundary condition poses a rele- 
vant perturbation to the original equation. While it is 
not possible to capture its effect by the simple methods 
mentioned above, it can be understood using physical ar- 
guments. The mechanism turns out to be very powerful 
and extends far beyond the EW problem. 

The Edwards- Wilkinson equation Q] describes the 
temporal evolution of an interface characterized by its 
height <p(x, t) over a substrate of length L, x € [0, L], at 
time t under the influence of a thermal noise T)(x, t). In 
one dimension it reads 

dt<t>( X) t) =Dd*4>(x,t)+r}(x,t) , (1) 

with a diffusion constant or surface tension D. The initial 
conditions are usually |J, |5j chosen to be 4>(x, t = 0) = 
and periodic boundary conditions (PBC) are applied. 
The central property one is interested in is the roughness, 
defined as 

w 2 {t,L) = (^)-(^ 2 ) , (2) 

where A denotes the spatial average, A = J Q dxA{x) 
and ( ) is the ensemble average, averaging over all real- 
izations of the noise r\. In order to determine w 2 , the 



only property of r\ which enters is 

(r](x, t) V (x',t')) = T 2 8(x - x')S(t - t') , (3) 

where T parameterizes the strength of noise. To fully 
specify the noise, the higher order correlations are usu- 
ally chosen to be those of Gaussian white noise and the 
average is set to Mx, t)) = 0. Assuming a Family- Vicsek 
scaling behavior |fj, three exponents, a, (3 and z, are de- 
fined for the asymptotic behavior of w 2 

w i i t,L)=aL 2 «g(^ , (4) 

with appropriate, system dependent parameters a and 
b which make the universal scaling function G(x) a di- 
mensionless function of a dimensionless argument. This 
function is supposed to behave like Q(x) cx x 2<3 for small 
arguments x and to converge to a non-zero value for large 
arguments. In the initial growth phase, w 2 is supposedly 
independent of L, so that (3 = a/z. 

Because of the small number of independent parame- 
ters, the exponents characterizing Q can be determined 
immediately, simply by using dimensional analysis. The 
roughness based on the ensemble of solutions of Q can 
be written as 

w 2 (t,L-D,T) = ^g^ , (5) 

with all independent parameters being listed on the left 
hand side. Eq. JSJ is the unique way of writing w 2 in 
the form prescribed by with all i-dependence being 
absorbed into a dimensionless function of a dimensionless 
argument. By assuming the existence of the appropriate 
limits or by simply comparing J^Jl to (@J), the exponents 
are determined to be: 

a = 1/2 (3=1/ A z = 2 (standard EW) (6) 

Of course, this can easily be confirmed by exact solutions 
00. 

The exponents JHJ remain unchanged, if fixed bound- 
ary conditions (FBC), (j)(x = 0,t) = <j>(x = L,t) = 0, 



2 



are introduced; the exact solution changes, but as the 
new boundary condition does not contain any new non- 
zero parameter and therefore cannot introduce any new 
scale, the scaling form Eq. JSJ) must necessarily remain 
unchanged; there are just not enough (independent) pa- 
rameters for the form Eq. (0 to change. 

On the face of it, the EW equation becomes much more 
complicated, if a drift or convection term is introduced 
into JU, 



dt<j){x, t) = Dd 2 4>(x, t) + vd x (j){x, t) + T)(x, t) 



(7) 



where v denotes the drifting velocity. Such a term is 
probably present in most experimental setups, for exam- 
ple due to the presence of a gravitational field and a small 
tilt of the substrate but it is not obvious whether the re- 
sulting drift is strong enough to have a significant effect, 
see also the crossover length calculated below). Here, it 
is worth mentioning the general role of a drift term as 
discussed in the context of anomalous aging by Luck and 
Mehta 0. 

In case of PBC, the solution of Eq. JJJ is related to the 
solution of Eq. Q by a simple Galilean transformation, 
4>'(x,t) = 4>{x + vt,t), which leaves the noise correla- 
tor Q unchanged. Because of translational invariance, 
the correlator (4>(x, t)<j>{x' , if)) (spatially) depends only 
on x — x' + v(t — if) and becomes independent of v for 
equal times t' = t. Because w 2 in J5J is a functional 
of the equal-time correlator, it is also independent of v. 
Thus, in case of PBC, the exponents remain those listed 
in Eq. ©, 

a = 1/2 = 1/ '4 z = 2 (EW with PBC and drift) 

(8) 

Dimensional analysis, however, leads to a scaling form 



w (t,L;D,T,v) = -=-0 



T 2 L „ tD vL 



D \ L 2 ' D 



(9) 



the behavior of which is unknown; as is shown below, 
Q(x,y) might behave like y" 1 / 2 for large x. 

It is very tempting and in fact leads to the expo- 
nents previously reported in the literature Q, to ar- 
gue that the drift term in (JJJ is relevant compared to 
the diffusion term. In fact, this is what an ansatz 
4>(bx,b z t) = b a (f>(xA) suggests. The idea behind this 
much-used method |9|,ll0|, is to determine the terms dom- 
inating the large scale behavior of Eq. J7J). It is often 
applied to determine relative relevancy of higher order 
terms in Langevin equations. In fact, the exact relation 

<j)(x,t;bL,D,T,v) = b^^x/b, t/b z \ L, Db z ~ 2 ,T, vb*^ 1 ) 

(10) 

seems to indicate that </> behaves on large scales like cj) 
on small scales with a diffusion constant D reduced by 



a factor b~ x compared to the reduction of v. Thus, the 
diffusion constant is expected to drop out for large system 
sizes. Repeating the dimensional analysis without the 
diffusion term, i.e. setting D = 0, then yields 

a = = z = 1 (EWd from coarse graining) . 

(11) 

This is consistent with a scaling law which can be derived 
from It states that any a derivable from dimen- 

sional analysis must obey Q 



a 



(12) 



Below it is shown that this is not the case. There is 
no reason to assume that <j) becomes asymptotically self- 
affinc. 

As shown above, this result, suggesting a smooth or 
flat interface || , is wrong at least in the presence of peri- 
odic boundary conditions, because the velocity v simply 
disappears due to translational invariance. The key ques- 
tion is then: Are the exponents Ijlll) recovered if the drift 
term cannot be gauged away? Does the coarse graining 
and dimensional analysis arguments presented above fail 
only because of translational invariance, which effectively 
removes v as free parameter from the problem? The an- 
swer to both of these questions turns out to be negative. 
The reason for that is an extremely efficient mechanism, 
effectively "wiping out" the stationary roughness. 

In order to answer the questions raised above, fixed 
boundaries are applied again, so that the drift cannot be 
gauged away. It is worth stressing that fixed boundaries 
correspond to a finite substrate or a finite region of ex- 
posure to the noise, i.e. these boundaries are much more 
natural than periodic boundary conditions. 

Eq. (0 with FBC can be solved using a saddle point 
approximation. In the following only a few technical de- 
tails of this calculation are presented. After expressing 
in a dimcnsionless form, the solution can be written 

as 



^(y, T ;<7)=/ d y' dT'<p (y-y',T-T')£(y',T') , 
Jo Jo 

(13) 

where y = x/L, r = t/(L 2 /D), q = vL/D and 
ip(y,r;q) = l/(Ty / L/D)(f)(x 1 1; L, D, T, v) as well as 
£(y,r) = l?l 2 /(T\^D)r](x,t) are dimensionless quanti- 
ties. The propagator tpo of this problem essentially con- 
sists of two Gaussians "wrapped around a circle of radius 
2" from a mirror charge trick, which are closely related 
to Jacobi's tf 3 -function IlllllJ. This sum is multiplied 



by an exponential e 2 
the drift: 



i Tq , which accounts for 
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(y-y +2n)* (y + y Q +2 n )' 

6 4r — - g 4t 



2 sin(n7ry) sin(n7r?/o) 



(mr) t e -%{y-yo)q--gTq 



(14a) 
(14b) 



n=l 



where by convention T2 > ti and y.Q is the initial position. The saddle-point approximation is required when the 
spatial averages © are taken. One finds 



w 2 {t,L) = T 2 x < 



t t \v\t^ 2 t 2 \v\ c ^ L 

, H for t < A 

2nD L 3LV2^D 2L 2 W 



1 



3^ V 2 I^I 



for i > ^ 



(15a) 
(15b) 



which becomes exact for divergent vL/D, i.e. especially 
in the thermodynamic limit. Inspecting the leading terms 
immediately gives 

a = 1/4 /3 = l/4 z = l (EW with FBC and drift) 

(16) 

This result is surprising, because these exponents are 
anomalous, as they do not correspond to what seems to 
be suggested by dimensional analysis: Not only do they 
contradict (JSJ, JSJl and (|11|) . they do not even obey (|12fl . 
In fact, neither the drift term, nor the diffusion term be- 
come irrelevant ; they are both crucial for (|16|l . More- 
over, the set of exponents switches from (jHJ to l|16|l if the 
boundary condition are changed from periodic to fixed 
in the presence of a drift term. The Edwards- Wilkinson 
equation with drift term and fixed boundary conditions 
is a linear problem displaying anomalous exponents. 

The fact that the above result cannot be reconciled 
with dimensional analysis might suggest that the mecha- 
nism leading to these exponents is very subtle. However, 
it turns out that it can be understood quite easily. 

The drift term makes the entire interface configuration 
move from one boundary to the other. Without noise, a 
peak starting somewhere in the bulk gets slowly moved 
by the drift to one of the boundaries, while diffusively 
broadening. It eventually disappears at the boundary. 
The time it spends between the boundaries depends on 
the starting position and the direction of the movement. 
The maximum time is L/v, which is also the maximum 
time, any noise-generated structure has to develop. How- 
ever, as known from the model without drift, it takes time 
L 2 /D in order to fully develop the roughness. Thus, for 
L 2 /D ^> L/v, i.e. L ^S> D/v = Lx, the interface will re- 
main in its initial growth phase; the characteristic length 
scale Lx represents an effective cutoff for correlations. At 
the same time it enables the system to display anoma- 
lous exponents. Depending on the direction of the drift, 
the interface "comes out" of the left boundary initialized 
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FIG. 1: A qualitative picture of an interface snapshot with its 
"local roughness", obtained in a numerical simulation; scales 
are irrelevant. Upper panel: An example of an interface pro- 
file with fixed boundaries and drift term. Lower panel: The 
ensemble averaged local width squared (numerical data, cir- 
cles) is proportional to the square root of the position where 
measured, x 1 ' 2 (fitted, dashed line). 



to <p = and moves to the right boundary, as shown in 
Fig- n The average "age" is proportional to tx = L/v, 
so that according to (QJ w 2 oc ti' 2 oc L 1 ! 2 , therefore 
a = 1/4. Regarding j3, the interface cannot "see" the 
drift initially, so that f3 — 1/4 just like for the case with- 
out drift, JBJ. Indeed, even the amplitude of the lead- 
ing term in (|15a|l corresponds to the amplitude obtained 
for the problem with periodic boundaries. The identity 
a = f3 already indicates z = 1, which can also be derived 
from the fact that saturation should be reached as soon 
as the interface has swiped through the system once, i.e. 
after tx — L/v. 

To test the validity of the results, Eq. I|15f) has been 
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compared to numerical simulations, based on a straight- 
forward integration of (JJJ using Euler's method. For sys- 
tem sizes not too small (L > 128), even the higher order 
corrections were reproduced. The mechanism is illus- 
trated in Fig. [2 The upper panel shows a snapshot of an 
interface configuration. The "local age" of the interface 
can be read off the local roughness (in an appropriate 
ad hoc definition) as shown in the panel below, because 
spending more time between the boundaries increases the 
local roughness according to w 2 oc t 1 ! 2 = (x/v) 1 ^ 2 , with 
x being the position where the roughness is measured. 

The physical explanation presented above goes beyond 
the EW equation; provided that the crossover time of 
the original model without drift scales faster in L than 
ix, i.e. z > 1, the argument should apply, so that at 
sufficiently large system sizes a obtains the value of (3 and 
therefore z = 1. It is a very efficient mechanism, which 
works under very general circumstances even in the most 
simple, linear case. It therefore speaks a clear warning 
as to the interpretation of numerical and experimental 
studies: the true value of a might have been "washed 
away" by a very small drift. 

The KPZ HllllEl equation (z = 3/2) is particularly 
interesting: Its non-linearity is only important during the 
initial growth phase [fi — 1/3) and becomes insignificant 
in the stationary regime (a = 1/2). However, with an 
additional drift term the equation should remain in the 
initial growth phase, with the non-linearity present in 
the stationary state. Indeed, using a Cole-Hopf transfor- 
mation the problem can be reduced to an equation 
similar to Eq. J7J). However, preliminary numerical tests 
did not fully confirm this correspondence. 

It is not yet completely clear how to generalize this 
arguments to higher dimensions. For two dimensions it 
is tempting to speculate whether exponents observed in 
experimental molecular beam epitaxy are related to such 
a drift term, for example when a « j3 |l4| or when a 
is close to /?kpz ~ 0.24 |lfj,|l5|. Interestingly, only one 
boundary needs to be fixed in order to observe the phe- 
nomenon, namely the boundary the velocity points away 
from, in Fig. ^ the left boundary. 

The mechanism stresses once more the relevance of 
boundary conditions as prominently pointed out by Lan- 
dau and Binder [lfij . However, it is worth emphasizing 
that in the present case, the change of boundary condi- 
tions leads to a change of the bulk critical exponents. 

Even though the physical explanation provides a very 
clear picture of the mechanism, it is not obvious how a 
drift term exactly affects other models, with, for example, 
quenched or conserved noise, with an additional term 
— rn<j) or with additional scaling laws. 

In conclusion we have presented a remarkably simple 
mechanism which reduces the roughness exponent to the 
value of the growth exponent for any small amount of 
drift in the Langevin equation in the presence of fixed 
boundary conditions, provided that in the original model, 



the dynamical exponent z is larger than unity. On suf- 
ficiently large scale, this mechanism should be visible in 
many experimental and numerical systems. Most unex- 
pectedly, it can even be found in the Edwards- Wilkinson 
equation, which consequently shows anomalous expo- 
nents, depending on the boundary conditions imposed. 
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